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In this paper we discuss two methods of resumming the leading and next to leading order in 
1/N diagrams for the quartic O(N) model. These two approaches have the property that they 
preserve both boundedness and positivity for expectation values of operators in our numerical 
simulations. These approximations can be understood either in terms of a truncation to the in- 
finitely coupled Schwinger-Dyson hierarchy of equations, or by choosing a particular two-particle 
irreducible vacuum energy graph in the effective action of the Cornwall-Jackiw-Tomboulis for- 
malism. We confine our discussion to the case of quantum mechanics where the Lagrangian is 
L(x,x) = (1/2) X^i=i ®i ~ (p/8JV) [X^ili 2 -? ~ r o] 2 - The key to these approximations is to treat 
both the x propagator and the x 2 propagator on similar footing which leads to a theory whose 
graphs have the same topology as QED with the x 2 propagator playing the role of the photon. The 
bare vertex approximation is obtained by replacing the exact vertex function by the bare one in the 
exact Schwinger-Dyson equations for the one and two point functions. The second approximation, 
which we call the dynamic Debye screening approximation, makes the further approximation of 
replacing the exact x 2 propagator by its value at leading order in the 1/N expansion. These two 
approximations are compared with exact numerical simulations for the quantum roll problem. The 
bare vertex approximation captures the physics at large and modest N better than the dynamic 
Debye screening approximation. 

PACS numbers: 11.15.Pg,11.30.Qc, 25.75.-q, 3.65.-W 



I. INTRODUCTION 

The need to understand quantum systems in real time 
in a quantum field theoretic setting arose from attempts 
to understand various early universe scenarios. These 
scenarios are based on the evolution of scalar fields ei- 
ther through their role as inflation fields or as topologi- 
cal defect forming fields. One would like to understand 
the quantum evolution of these fields rather than rely on 
unjustified treatments based on studying their classical 
evolution. The study of the "slow rollover" transition in 
an upside down harmonic approximation by Guth and 
PiQ] was the first attempt to understand whether clas- 
sical approximations could be justified. However, one 
really needed to go beyond the harmonic approximation 
to address the nonlinear aspects of double well (and Mex- 
ican hat) potentials. These non-linear aspects effect pro- 
duction of topological defects as well as the nature of the 
oscillation at the bottom of the well which causes reheat- 
ing. 

Our ultimate goal is to be able to describe accurately 
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over relevant time periods the nonlinear aspects of quan- 
tum field theory evolutions. Although in one-dimensional 
quantum mechanics, one can rely on a numerical solu- 
tion of the Schrodinger equation to understand the time 
evolution of the system accurately over long time peri- 
ods, in field theory contexts the numerical solution of the 
functional Schrodinger equation is presently beyond the 
reach of the largest computers. One important question 
is how to decrease the number of degrees of freedom in 
a manner consistent with certain physical requirements 
such as conservation of energy, preservation of positiv- 
ity and boundedness of expectation values. Although 
this is guaranteed in variational approximations, approx- 
imations based on various truncation schemes, whether 
perturbative or non-perturbative in nature often fail to 
preserve these physical requirements. For example, naive 
truncations of the coupled Green functions equations be- 
yond the truncation at the two-point function level lead 
to secular behavior (unboundedness at late times). This 
is also true for the 1 /N expansion which is derivable from 
an effective action. The second question is, after guaran- 
teeing these properties, how accurately have we described 
the time evolution. 

The simplest truncations of the field theory have been 
based on gaussian variational methods^, ||, or the re- 
lated leading order in large-iV approximation (LOLN) 
|Q, 0]. These two approximations can be shown to be 
equivalent to a classical Hamiltonian dynamics for the 
variational parameters (or equivalently the Green func- 
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tions) which leads to probability conservation at the 
quantum level so that the results always lead to con- 
served energy, and positive and bounded expectation val- 
ues. Unfortunately, hard scatterings which lead to ther- 
malization are ignored so that important physics is left 
out. The approximation also is numerically inaccurate 
after a few oscillations in quantum mechanical applica- 
tions, unless the anharmonic coupling constant g divided 
by the number of fields N in an 0{N) model is quite 
small. In this paper we will be comparing our methods 
of going beyond mean field theory (Hartree or large- N) 
with exact numerical simulations of a quantum mechani- 
cal O(N) model. In this way we can see how accurate the 
approximations are as a function of N as well as study 
numerically if the approximation maintains the various 
physical requirements we posit, such as boundedness and 
positive definiteness of expectation values. The reason 
for using this quantum mechanical model is that exact 
simulations can be done at all N, so that accuracy of 
the method as a function of the parameter 1/N can be 
studied. By restricting ourselves to a quantum mechan- 
ics problem we unfortunately will not be able to study 
questions of thermalization. A complementary approach 
has been undertaken by Aarts, Bonini, and Wetterich|| 
where they consider classical 1+1 dimensional 4> 4 field 
theory (for N = 1). There one can look at some aspects 
of classical thermalization (as long as one keeps a cut- 
off because of the Raleigh- Jeans divergence) but one is 
restricted to low values of N so one cannot study the 
N dependence of the result. Also one cannot study the 
quantum aspects of the problem. In the above paper, 
Aarts et. al. study a truncation of the Green functions at 
the four-point level, which is known to lead to unbound- 
edness and secularity in quantum mechanical (as well as 
classical) applications. It will be interesting in the future 
to apply the approximations we are using here to clas- 
sical 1+1 dimensional 4 to see if, and how well, they 
describe the thermalization. 

There are several ways of approaching the problem of 
thinning the degrees of freedom of the quantum field the- 
ory. One of the earliest was based on making a variational 
approximation to the functional Schrodinger equation. 
The variational approach has the advantage of leading 
to a Hamiltonian dynamical system for the variational 
parameters as well as to a density matrix which has pos- 
itivity properties. Energy conservation and positivity 
and boundedness of expectation values are automatically 
guaranteed. However, even for the simple problem of the 
quantum roll, the gaussian, or time dependent Hartree 
approximation, studied by Cooper, Pi and Standoff , 
and improvements which are based on trial wave func- 
tions of the form of a polynomial times a gaussian[Q, 
were found to be only accurate for relatively short time 
periods (one or a few oscillations) when compared to the 
exact numerical solution of the Schrodinger equation. In 
quantum mechanics, except for exceptional situations, 
the wave function in multiwell situations gets very com- 
plicated very quickly and is not easily described by a 



small number of variational parameters. 

A second approach has been a direct 1/N expansion 
of the path integral in the Schwinger-Keldysh-Bakshi- 
Mahanthappa closed time path formalism[||. In this ap- 
proach the connected Green functions have the property 
that they start at order G211 oc l/N'' 1 ^ 1 . Thus if we 
retain only a certain order in the expansion, there is 
a truncation in the order of Green functions retained. 
This approach was applied recently to the quantum roll 
problem ||] and was found to suffer from the secularity 
problem — although the short time behavior of the re- 
sult was improved by including 1/N corrections, an exact 
reexpansion in terms of 1/N leads to corrections in the 
Green functions that are of the form ±t/N and so the 
individual corrections become unbounded as well as non 
positive definite. In this approach, although energy is 
conserved, individual contributions are not positive defi- 
nite and unphysical behavior is found. 

A third approach has been to consider the complete 
set of equal time Green functions. These obey first order 
local equations in time, as in the Schrodinger approach. 
This approach has been nicely systematized and an equa- 
tion for the generating functional obtained by Wetterich 
and collaborators in a series of papers |l(J. However, 
naive truncations of the equal time Green function hier- 
archy again have the problem that although there is a 
conserved energy, one cannot show that this truncation 
(except at the two-point level) corresponds to a posi- 
tive definite probability so that expectation values are 
not necessarily bounded or positive definite. Truncated 
at the two point function level, this approach is identi- 
cal to the Hartree approximation. However, simulations 
based on truncations assuming 6th order or 8th order 
1-PI graphs, could be set to zero, were carried out for 
the O(N) x 4 th oscillator problem, and secularity was 
found for many choices of initial conditions |]llf . So we, 
as quantum field theorists, having entered the domain 
of nonequilibrium phenomena, are now beset with all 
the problems faced by our plasma and condensed matter 
brethren more than 40 years ago! 

In both quantum and classical many-body systems, the 
dynamical equations are an infinite hierarchy of coupled 
equations which relate given ensemble averages, whether 
classical or quantum, to successively more complicated 
ones. To make the solution of this hierarchy possible, 
some truncation scheme is necessary. Most naive trun- 
cation schemes which, for example, just truncate the hi- 
erarchy of coupled correlators at a particular order, do 
not preserve various physical properties required of the 
system — such as positivity of the spectral components 
of the Green function and conservation of probability. A 
corollary of this is that in most perturbation schemes, 
secularity arises quickly with each term in the pertur- 
bation series, growing with higher powers of the time 
t. In his seminal paper of 1961, Robert Kraichnan[fl2| 
discussed in detail the key issues and obtained a partial 
solution to the problem by demanding that the approxi- 
mations one should use should correspond to some physi- 
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cally realizable dynamical system. This would guarantee 
positivity and secularity would be avoided. The reason 
why variational approximations avoid these problems is 
exactly because they lead to a Hamiltonian dynamical 
system for the variational parameters (which are related 
to equal time correlation functions). He also discussed 
scenarios where particular classes of graphs, which con- 
tained the relevant dynamics, are summed and he sug- 
gested some physically motivated approximations which 
did not suffer from any diseases. In field theory one rarely 
has the parameter control to make such guesses, how- 
ever some progress in QCD has been made by summing 
hard thermal loops |l3| , which already tells us some of the 
graphs that we want to include. In plasma physics, one 
wants to make sure that the approximation to the dy- 
namics is robust enough so that the photon propagator 
includes polarization effects, which give Debye screening. 
This is related to the hard thermal loop summation in 
QCD. 

To find resummation schemes that avoid the secularity 
problem we will rely on the experience of our many-body 
and plasma physics friends. To calculate the conductivity 
of a non-relativistic plasma, it is known what graphs are 
necessary to sum in order to get agreement with experi- 
mental results^], . Basically the conductivity is found 
from the vertex function which must satisfy an integral 
equation which sums ladders of the Debye screened pho- 
ton propagator. The two approximations we will discuss 
here will differ on whether the equivalent of the Debye 
screened photon propagator for the anharmonic oscilla- 
tor is treated in lowest order in mean field theory, or is 
self-consistently determined. In studying the conductiv- 
ity of a relativistic plasma the first approximation has the 
advantage of obeying the correct Ward identities (but vi- 
olating energy conservation to order 1/N) whereas the 
second preserves energy conservation but violates Ward 
identities (to order 1/N 2 ). Here we are not studying 
QED, and the Ward identities of the 0(N) model for the 
quantum mechanics problem are much simpler than those 
of QED and energy conservation is a more important con- 
straint on the accuracy of the answer. We will include 
both approximations here mainly because of the recent 
interest in the gauge invariant approximation for the rel- 
ativistic plasma p9|, and also because in truncations of 
Schwinger-Dyson equations, it is often too difficult to 
solve for the photon propagator self consistently, and so 
one is often forced to try the more drastic approximation 
of using the mean field propagator in the resummation 
scheme. By studying this approximation in a quantum 
mechanics problem we will see the shortcomings of such 
an approach. 

In what follows we will discuss two approaches to ob- 
taining the above two truncations of the exact Schwinger- 
Dyson equation and apply them to the problem of the 
quantum roll — the long time behavior of N coupled 
anharmonic oscillators with "radial" symmetry in an N- 
dimensional space. This particular problem has been 
studied by us previously M exactly and in the next to 



leading order in the large-iV approximation (NLOLN) 
and is interesting because exact numerical solutions can 
be found for arbitrary N. What we found previously, is 
that for the parameter set studied (g ~ 1, M 2 = 2), the 
next to leading order in large-iV contributions became 
unbounded for N < 21. For larger N, where the ap- 
proximation was physical, it had the failing that it was 
unable to track the spreading of the exact wave function 
which led to the envelope of the oscillations found for 
(x 2 (t)) contracting at late times and then reexpanding. 
A related study of large- N for quantum mechanics in the 
context of the equal time correlators by Bettencourt and 
Wetterich|nj, also displayed growing modes for various 
initial conditions. 

The resummation presented here will allow one to track 
the contraction for some period, but at later times it also 
fails in that it leads to small oscillations about a fixed 
point value. In field theory settings, where one hopes 
that this approximation will lead to thermalization, op- 
timistically this fixed point behavior will become physical 
and be related to thermal equilibration. Whether this is 
true or not can be checked by studying this approxima- 
tion for classical evolutions averaged over a distribution 
of initial conditions described by a an initial probability 
distribution in phase space. 

In what follows we will present numerical solutions for 
the quantum roll problem for the O(N) model, and com- 
pare them to these two different approximations to the 
Schwinger-Dyson equations, which sum infinite numbers 
of leading order and next to leading order in 1/N graphs. 
Our approach will be to introduce a composite "field" 

X = | a? - r^j , 

which is treated on equal footing to the field x. By doing 
that, the Schwinger-Dyson equations for the theory will 
have the same topology as those of QED with x play- 
ing the role of the electron and \ the role of the photon. 
At leading order in large N in N- flavor QED, one sums 
all the fermion loop vacuum polarization corrections to 
the photon propagator which gives the Debye screening. 
Here the bare photon propagator is replaced by a local 
interaction in the graphs for the \ propagator in LOLN. 
The next consideration, important for charged plasmas, 
is that to obtain reasonable agreement with experiments 
on the conductivity of the plasma, the vertex function 
must sum all the ladders with the Debye screened propa- 
gator as the kernel in the integral equation. The two re- 
summation schemes which we discuss in this paper both 
have this property. 

The approximation which we call the bare vertex ap- 
proximation (BVA), uses the full Green function for x as 
well as the full Green function for % in a 2-PI Hartree 
graph contribution to the effective action. This is in con- 
trast to an earlier scheme for going beyond 1/A|l4| using 
the 2-PI formalism which is based only on the x Green 
functions. The BVA approximation sums an infinite Ge- 
ometric series of 2-PI graphs of the single field formalism. 
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Recent simulations in a toy 1+1 dimensional scalar field 
theory|l5) show that the approximation described infill 
already has the ability to thermalize arbitrary initial con- 
ditions, so we are confident that the BVA approximation 
will also have that feature when applied to a field the- 
ory problem. The BVA can also be obtained by setting 
the full vertex function to unity in the Schwinger-Dyson 
equations for the one- and two-point functions with exter- 
nal sources hence the origin of its name. The second ap- 
proximation we will study, which we call the dynamic De- 
bye screening approximation (DDSA) , makes the further 
assumption that the full x propagator can be replaced 
by the lowest order in 1/N composite field propagator 
in all the integral equations. The main interest in the 
DDSA results from it being the lowest order resumma- 
tion scheme that exactly preserves QED Ward identities. 
Both these approximations are free from the difficulties 
found in the perturbative l/N expansion, which we dis- 
play for comparison. We find that the BVA is accurate 
at modest times < 25 oscillations when N > 10. At later 
times it settles down to oscillating about an unphysical 
fixed point. The DDSA approximation violates energy 
conservation at order 1/N and as a result becomes inac- 
curate after several oscillations. In spite of this, it is nu- 
merically more accurate for a longer period of time than 
the Hartree approximation at small and modest values of 
N. 

It should be kept in mind that quantum mechanics and 
quantum field theory are very different. For example, in 
the quantum mechanics application discussed here, the 
graphs of the 0(1/N) corrections do not correspond to 
interparticle collisions (as they do in field theory) since we 
are restricting ourselves to one-particle quantum mechan- 
ics. Nevertheless quantum mechanical examples provide 
excellent test beds for key issues such as positivity vio- 
lation, boundedness, and late time accuracy of the ap- 
proximations. It is precisely these questions that we are 
hoping to understand in this paper. 



II. THE O(N) MODEL 

The classical Lagrangian for the O(N) model of N non- 
linear oscillators is given by: 



JV 
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The Schrodinger equation for this problem is given by 

dijj{x, t) 



1 N a 2 1 



2^ dx 

1=1 



dt 



(2) 



where V(r) is a potential of the form 



V(r) = — 
V ' 8N 



9 /„2 2\ 2 2 



'07 > 



N 

= £*?• (3) 



For the quantum roll problem there is spherical symme- 
try. This means that we can assume a solution of the 
form ip(r,t) = <p(r,t)/r( N ~ 1 " 2 , in which case the time 
dependent Schrodinger equation for <fr(r, t) reduces tojjyj): 



29^ 



+ U{r) \ 0(r, t) = i 



d<t>(r,t) 

at 



(4) 



with an effective one dimensional potential U(r) given by 
(W-l)(jV-3) , g 



U(r) 



y i i 2v 

SAT V 0/ 



(5) 



sr" 8N 

It is this equation that we will solve numerically to obtain 
exact numerical solutions as a function of N. U(r) has a 
minimum at r — r m \ n . In our simulations, we have fixed 
our mass scale M 2 , defined as the second derivative of 
U(r) at the minimum, to have a value of 2, independent 
of V. 

Returning to the Lagrangian formulation, it is useful 
for the purposes of obtaining a large- ./V expansion to in- 
troduce scaled variables: 

(6) 



'Nxi , r -» V Nr . 
Then the Lagrangian scales by a factor of N: 



L/N = L N (x,x) 



1 



i=i 
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N 



(7) 



We use these scaled variables in this paper, so that the 
rescaled ro ~ 1. Next we introduce a composite coordi- 
nate x by adding to (0) a term: 



1 

29 



f(x>- 

1=1 



(8) 



The Lagrangian (|7|) then becomes: 
1 



L N (x,x;x,x) = ^2 



roX 



+ ^ + Jx, (9) 



where we have also added sources ji and J coupling to 
Xi and x respectively. From this Lagrangian we get the 
Heisenberg equations of motion for the operators Xi(t) 
and x(t): 




X(t)xi{t)=ji(t), 



Y,xKt)-r 2 o)-J(t) 



(10) 



Here, and in the following, we indicate operators by 
"hats." Taking expectation values with respect to an 
initial density matrix we obtain the c- number equations: 



(*(*)> 



(£<(*)) + (*(*)&(*)) 

N \ 



£*?(*) 




(ii) 
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By rewriting the quartic interaction in terms of the com- 
posite field x, the induced interaction of the form \x 2 is 
reminiscent of N flavor QED with interaction A^ij^tpi. 
The fact that these two theories have the same topologi- 
cal structure will allow us to use the intuition gained in 
classical plasmas to make appropriate approximations. 

To simplify notation we include all independent coor- 
dinates in one vector. We define: 



x a (t) = [x(t),x 1 (t),x 2 (t), . . . ,x N (t)] 
j a {t) = [J{t),ji(t),j 2 (t),... ,j N {t)]. 



(12) 



for a = 0, 1, . . . ,N, and where J(t) = J(t) - r%/2. Ab- 
sorbing the factor r§/2 into the current means that J(i) 
is not zero when J(t) is set to zero. Greek indices run 
from to N, whereas Latin indices go from 1 to N. Us- 
ing this extended notation, the generating functional Z[j] 
and connected generator W[j] is given by the path inte- 
gral: 

N , 

z[j] = e iNW[j] = Y[ dx a exp{iNS N [x;j]] (13) 
where the action Sn[x;J] is given by: 

S N [x;j] = ~Y / [ dt [ dt'x a (t)A-^[x}{t,t')x (t') 

a,p C ^ C 

+ J2 I dtx a (t)j a (t) , (14) 

and where A~^[x] (t, t') is given by: 

* l r i/ /\ fD- l (t.t') \ 

A^N(t,o = ( i' ctfwwy ' (15) 



with 



1 



D- 1 (t,t') = —S c (t,t') 



9 

G^[ X }(t,t') = {^ + Mc(*,0 . (16) 

In what follows it will be useful to introduce another 
matrix inverse Green function G~p[x](t, t') as follows: 



S 2 S N [x;j] 



(17) 



Sx a (t) 8xp(t') 

'D-\t,t') Kj\t,t')\ 
K K7\t,t>) Grj(t,t')J ' 

with D -1 (t,t') and Grj(t,t') given by Eq. @, and 
K- 1 [x){t,t') = K- 1 [x]{t ) t')=x l {t)5 c {t,t'). 

III. THE SCHWINGER-DYSON EQUATIONS 

The Schwinger-Dyson equations are integral equations 
for the Green functions. The Green functions can be 



obtained by functional differentiation of the path inte- 
gral for the generating function in the presence of exter- 
nal sources. After setting the external sources to zero, 
one obtains an infinitely coupled hierarchy of coupled 
equations for the Green functions. For an initial value 
problem, the boundary conditions on the Green func- 
tions can be implemented by using a time ordered prod- 
uct where the time ordering refers to the closed time path 
contour of the Schwinger-Keldysh-Bakshi-Mahanthappa 
formalism||. A detailed discussion of that formalism as 
applied to implementing the 1 /N expansion for this par- 
ticular problem is described in ref . || . One way to gen- 
erate the equations is to consider the identity fl7j]: 



n 



d.r 



Sx a (t) 



JNS N [x;j] _ 



0, 



from which we find: 



(18) 



J(t), 



— + X (t) f Xi(t) + -)d(t, t)/i = jt(t) , (19) 



where Xi{t) and x(t) are average values of the operators, 

-w-^lf -*<*»■ 
, (t)s ™ -<«.», 

and where the Green functions Q a p\j](t, t') are defined 
by: 



S 2 W[j] 



_ 5x a (t) _ 

" 6j(,(f) ~ 5j a (t)8j (t>) 

_ (V{t,f) Ki(t,t)\ , , 

~ (iCi(t,t') GiAt,?)) ■ [ZV) 



GoMM 



Eq. ( |L9| ) is identical to Eq. ([Ll|). In this equation and in 
what follows, Xi and x now correspond to the expectation 
values: 

The Green functions are explicitly given by 



V(t,t>) = 

iCi{t,t') = 



S 2 W[J,j] 
SJ{t) SJ{t') 

S 2 W[J,j] 
Sji(t)8J(t') 



TCi{t,t') = 



Gi,j(t,f) 



S 2 W[J,j] 
S 2 W[J,j] 



The integrability conditions require that K-i(t,t') = 
JCi(t',t). To obtain the Schwinger-Dyson equations it is 
advantageous to Legendre transform to the expectation 
value of the coordinate variables x a (t), as the indepen- 
dent variable instead of the currents. The effective action 
generating functional of 1-PI graphs is given by a Legen- 
dre transformation: 



t[x] = w[j] - f dtJ2Mt)Ut)} 

Jc 



(21) 
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So since j a (t) — —Sr[x]/6x a (t), the equations of motion 
( |l9| ) give values for derivatives of T [x] : 

8T[x] 1 , , 
0X{t) 9 



He 



ST[x] 

Sxi(t) 



x 2 % {t) + -gu{t,t)/z -rl\ (22) 



^+ x (t)\x i (t) + ±!C i (t,t)/i. (23) 



However the Green functions here, Qn{t,t) and fCi(t,t) 
are defined in Eq. (|2(]) as functionals of the currents j a (t) . 
These must be expressed as functionals of x a (t) by in- 
verse relations. We define these inverse Green functions, 
which arc functionals of x a (t), by: 



5x a (t) 5x p {t') 



/C-\t,t') grUt,t') 



where explicitly 

T>-\t,t') = 
ici 1 (*,*') = 



s 2 r[ x , 



S X (t) 5 X (t') 



S 2 T[ X ,. 



6 X {t)6xi{t>) ' 
5 Xi {t)5 X {t') ' 



S 2 T[x,- 



6xi(t) Sxj(t') 



Again we have K i (t, i') = IC i (t', t). The inverse Green 
functions are given by differentiating the equations of mo- 



tion, Eqs. ( p2[ ) and (|23J), with respect to the coordinates. 
Using 



we find: 

&r 7 (i 3 ) 



di 4 / di 5 y^^ a ^[j](ti,i 4 ) 



X T Sten [x}(t4,t 5 ,t 3 )g e ^[j](t 5 ,t2) , 

(24) 

where V a ^ n [x\{t\, t2, t 3 ) is the three-point vertex func- 
tion, defined by: 



T a ,p, 1 [x](ti,t2, t 3 ) = 



5xj(t 3 ) 



S 3 T[a 



8x a (ti) Sxf3(t 2 ) Sx^(t 3 ) 
Explicitly, we find an equation of the form: 



(25) 



= G-^(t,0 + S a ^(t,0 , (26) 



where G^it^t') is given by Eq. (|l7|). The generalized 
self energy E a p{t,t') is given by: 

and where the polarization II(i, t'), self energy Ey(i, t'), 
and the off diagonal terms Sli(t, t') and f2j(t, f) are given 
by: 



n(t,t') = E / d< l/ d*2Si,a(*,tl)r a , j 9,o(tl,*a,* / )^,i(*2,t), 

Sy(t,t') = 4y] / d*l/ dt 2 ft,a(*,*x)r ajW (ti,*2,*')e/J,o(*2,*) 

n i(*>*') = ^y] / d*i/ d*2ei,a(t,ti)r a ,/>,i(ti,*a,t / )S/Sj(*a,*) 

a,/3,o(il,^2, t')Gpfi(t2, t) . (28) 



I 

In order to solve the equation for the two point function, hierarchy of coupled Green function equations that cor- 

Eq. (^6|) , one requires knowledge of the three point func- responds to solving exactly the Schrodinger equation, 
tion, defined by Eq. (^5|) . This in turn requires knowledge 

of the four-point function, ad infinitum. It is this infinite The matrix i nve rsion of Eq. (M) gives the set of cou- 
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pled equations, 



x 1 (t). We find: 



Q a ,pW) = G a ,f i ttJ)-Y\ f dtt [ dtaGa, 7 (t,ti) 



7,(5 

X%(tl,t2)Mt2/), (29) 



where 
with 

E 



Ga,p(t, t) — 



D(t,t') Ki(t,t') 
Ki{t,t') Gij(t,t') 



(30) 



d 2 



6ij+gxi(t) Xj (t)}G jk {t,t') 



(31) 



= S ik Sc{t,t') , 
D(t,t') = -gS c (t,t')+g 2 Y,^(t) G^t') Xj (t') , 

ij 

(32) 

Kj(t,t') = K,(t',t) =gY J G 3l {t,t')x t {t l ) . (33) 



When Xi(t) ^ 0, one notes that D(t, t 1 ) is not the inverse 
of t'). 

The vertex function Y a .j3^[x\{ti, t 2 As) defined in j25| ) 
is obtained by differentiation of Eq. (|2q) with respect to 



r a ,/3, 7 [x] (ti,t2,*3) 



fr) 

&E 7 (t 3 ) 



= fa,p,"t Sc{h,h) Se(ti,t 3 ) + ^ a ,p n [x]{ti,t 2 ,h) . (34) 

Here f i: j t0 = /o,;,j = / 4 ,oj = otherwise / is zero. 
&a,/3,~/[x]{ti ,t 2 ,t 3 ) is given by derivatives of the self- 
energy matrix: 



* r iu + + \ SE a [x](t 1 ,t 2 ) 



and is of order l/N. 

We are interested in resummation schemes that are 
exact to order l/N for (xf). We see from Eqs. ([}4]) and 
( pB| ) that it is consistent to replace r n g^[a:| (^1,^2^3) in 
Eq. (|2^) by the first term in Eq. (|3J) to obtain a re- 
summation which is exact to order l/N. To simplify 
our discussion of the exact Schwinger-Dyson equation for 
the vertex function, we will only consider the case of the 
quantum roll where Xi(t) = 0. 

Following the treatment of ref. |l8|], we have for the 
3-x vertex: 



A(ti,t 2 ,t 3 ) = — x ^\ t2) = V / dtj &hg lJ {t i M)G l k{hM)M 3k {t Al t 2 \h,t 1 ) , 

where Mj k (t4,t 2 ;t5,ti) is 1-PI in the channel x + x — > x + X- The lowest order in l/N contribution to M(ti 1 t 2 , t 3 ) 
is: 



M jk (t4,t 2 ;t 5 ,ti) = 5c{ti,t 2 )5c{h,ti) Qjk(h,ti) 
When Xi{t) — 0, the exact Schwinger-Dyson equation for the \~ x ~ x vertex is 



(36) 



r»_j (ti ,t 2 ,t 3 



SijSc(t 1 ,t 2 )5c(t 1 ,t 3 )- / dU dt 5 dt 6 dt 7 x 



^ r fe; (^4, t 5 ,t 3 ) Gkm{ti,h) Qin{h,t 7 ) K.i mn (t 5 ,t 2 ;t 7 , ti) + A(i 4 , t 5 , t 3 ) 2?(t 4 , i 5 ) V(t 6 ,t 7 ) IC 2l j(t 5 , t 2 ;t 7 ,ti) > . (37) 

klmn J 

I 



where ICi and ZC 2 are the s-channel 2-PI scattering ampli- 
tudes for the reactions: x + x — > x + x and x + x ~~ * x + x, 
respectively. 

This is shown pictorially in Fig. |l|. In general one then 
has to obtain equations for the 2-PI scattering ampli- 
tudes as well as for A. These will depend on even higher 
n-point functions, ad infinitum. In our approximations 
made at the two-point function level, the 2-PI s-channel 



scattering amplitudes K\ and K 2l used in the equations 
for the vertex function, will turn out to be graphs for 
one-particle exchange in the t-channel of the \~ an d x- 
particles respectively. 

In our truncations of the Schwinger-Dyson equations, 
we will always replace the full three-point vertex func- 
tion by the bare one in the equations for x and Q in 
the presence of external sources. Once this truncation 
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FIG. 1: Schwinger-Dyson equations for the vertex function. Solid lines represent the Qij(t,t') propagator and heavy wiggly 
lines are the T>(t, t') propagator. 



is made, then for the problem we are addressing here 
(the approximate time evolution of N quantum anhar- 
monic oscillators) one never needs any of the N point 
functions beyond the 1 and 2 point function equations. 
What will distinguish a further approximation we will 
call the DDSA is that we will also further approximate 
the x propagator to be that of the LOLN approximation. 

By making this bare vertex approximation in the equa- 
tions for the one- and two-point Green functions, we have 
not relinquished our ability to calculate in this approxi- 
mation all the higher connected Green functions. These 
are obtainable by further functional differentiation of the 
effective action. In particular if we wanted to use linear 
response theory (the Kubo formula) to obtain the electri- 
cal conductivity for a QED plasma, one would function- 
ally differentiate the equation for the inverse two-point 
for the electron function with respect to A^. In our prob- 
lem the photon is replaced by the composite field \, and 
the electron by Xi. 

Because of recent interest in studying plasma conduc- 
tivity in both QED and QCD, we will spend extra time on 
comparing the equations obtained for the vertex function 



in the three approximations considered here. In conduc- 
tivity calculations, it is necessary to sum all the ladder 
graphs in the equation for the vertex function to get good 
results for dilute plasmas. We will find that in NLOLN 
the vertex function is not an integral equation but is 
rather the sum of a few diagrams whereas the other two 
approximations lead to integral equations that sum an in- 
finite number of diagrams. Another issue is in preserving 
Ward identities. One of the reasons the large- N expan- 
sion was so interesting is that it is a complete reexpan- 
sion of the field theory which preserves Ward identities 
at each order. The QED plasma conductivity problem 
people Jl!| became interested in the DDSA because it 
exactly obeyed the Ward identities, whereas the BVA 
approximation violates Ward identities at order 1/N 2 . 
It is for this reasons we thought it appropriate to study 
the DDSA approximation, even though it violated energy 
conservation already at order 1/N, hoping that at least 
at large N it would be numerically accurate and satisfy 
Ward identities in QED applications. 

The exact formula for the energy is given by: 



E / N = \\Y. + *w w - to® - * 2 (*W 



(38) 



When Xi(t) = (xi(t)) = and ii(t) — 0, one obtains: 



E/N= 1 -^ d2Q ^ i 



dtdt' 



+ x (t)g(t,t)/i- rg X (t) - - 
t'=t 9 



x 2 (t) + -v(t,t)/i 



[ dtl [ dt2 [ d *3 2 , (*i,*2)0ii(*i,t)a«(t,ti)r ifc (t 1 ,t 3> ta)] > 

ijk J C ^ C J C > 



(39) 



where V jk{t\,tz,t<i) is the full vertex function given in 
Eq. @. 



IV. EFFECTIVE ACTION FOR 
TWO-PARTICLE IRREDUCIBLE GRAPHS 



Since the approximations we are going to consider have 
a simple interpretation in terms of keeping a particular 
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2-PI vacuum graph in the generating functional of the 
2-PI graphs, we would like to review this formalism fol- 
lowing the approach of Cornwall, Jackiw, and Tomboulis 
(CJT)@. 

The first Legendre transform of the generating func- 
tional W[j] of connected Green functions is widely known 
and used and is called the "effective action." The higher 
Legendre transforms (second, third, etc.) were intro- 
duced by De Dominicis and Marti n[f21[ in quantum statis- 
tics. Dahmen and Jona-Lasiniop2|7 and later Visil'ev 
and Kazanskiij2^], extended these ideas to quantum field 
theory. These methods were then used by Cornwall, 
Jackiw, and Tomboulis to discuss dynamical symmetry 
breaking in Hartree type approximations which later led 
to the second Legendre transformation formalism being 
called the CJT formalism. These higher order Legendre 
transformed actions have the advantage of being able to 
treat higher order Green functions on the same footing 
as the coordinates. 

We will first summarize the general results of that pa- 



per before proceeding to the specific approximations we 
consider in this paper. The method of CJT is to in- 
troduce one- and two-body sources for the coordinates 
x a (t) and the Green functions G a ,/3(t,t') in the action, 
and then make a Legendre transformation to the one- 
and two-point functions. The resulting action, as a func- 
tion of x and Q, contains a term which is the sum of all 
two-particle irreducible vacuum graphs. This term can 
be written using the vertices of the interaction and Q. 
We use the extended notation for the coordinates and 
one-body sources, given in Eq. (fHf). 

Thus, the generating functional Z[j, k) for the CJT ac- 
tion is given by: 

N 

Z\j,k] = e' w M = H dx a exp{iNS N [x;j,k]} , 
with 



where 



S N [x;j,k] = S class [x] + V / dtx a (t) j a (t) + ]- V / dt dt' x a (t) k a ,p(t, t') x0) , 

S c iassM = -i^ [ dt [ dt'x a (t)A-^[x}(t,t , )x l3 (t') = So + S in t[x}, 

/ dt [ dt'x a (t)A-^(t,t')x (t'), 

\ [ dtx(t)^2x1(t) , 
J c 



So 

Sint [x] 



(40) 

(41) 
(42) 
(43) 



and where A Q * Jt, t') is given by: 



£> -1 (t,i') 







G^t') = \ — \5 l Mt,t') 



effective action: 



T[x,G\=W\3,k\-Y i I dtx a (t)j a (t) 

\j2 j dt j dt ' Xf3(t')+Ga,f3 {t,t')} 



CJT showed that r[x,C/] can be obtained as a series ex- 
pansion in terms of 2-PI graphs. That is, introducing the 

with D~ l (t,t') given by Eq. (|l^). In this formalism, we functional operator, 

have separated out an "interaction" term, Eq. ( [i"3"|) , which 

depends on the coordinates x a {t), from a bare Green 

function Gg~^(t,t'), which is independent of the coordi- 
nates x a (t), in contrast to our previous definitions in 

Eq. (|l6|). The term rgx(i)/2 has been absorbed into the 

definition of the current J(t) in Eq. (p 



G~>](M') = 



S 2 S \x] 



5x a (t) 6x p (t>) 

D-\t,t') Kj^x^t') 
K^[ x ](t,t>) Grj[ x ](t,t')J ' 



(44) 



The second Legendre transform of W\j, k] is the CJT which is the same as the G a \[x](t, t') as defined in 
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FIG. 2: The vertex function for the BVA. The top figure represents Eq. (^) and the bottom figure represents Eq. (^8|). Solid 
lines represent the Gij{t,t') propagator and heavy wiggly lines are the T>(t,t') propagator. 



Eq. (|l7|), one can write the effective action in the form: 



in the exact Schwinger-Dyson equations for the self- 
energies, Eqs. (Pa). This gives for the BVA: 



T[x,G] = S c]aBS [x] + -Tr{ IniQ- 1 ]} 



'-Tr{G^[x}g -1} + T 2 [x,g}. (45) 



(47) 



The quantity ^[x, Q] has a simple graphical interpre- 
tation in terms of all the 2-PI vacuum graphs using ver- 
tices from the interaction term. The Hartrcc and leading 
order in large- N approximation for the x A potential was 
obtained by CJT using a single two-loop vacuum graph 
in the 0(N) theory written in terms of only the coordi- 
nates Xi. Our strategy for obtaining a resummation of 
the large- N approximation is to first rewrite the theory 
in terms of the composite field x, and the equivalent La- 
grangian given in Eq. (H). Using these new variables, we 
then choose for T^x, Q\ the 2-PI graphs shown in Fig. [|, 
which is now written in terms of the full \ an( i x propa- 
gators and the trilinear coupling x(i) xf(t)/2. 



i 

N 



s « (*> = jj (*» *i (*', *) + Sa (*, *)} . 



where we have used the symmetry property, &j(t, t') = 
Qji{lf, t) and Ki(t, t') = K.i(t\t). Thus we find fli(t, t'Y= 
fli(t', t). The self-energies (|47|) are then used in Eqs. ( J2q ) 
to find the one- and two-point functions. For the Green 
functions, wc find: 



£aj?(*> t ') = G a,p( t ' t ') + S BVA Q ,/3(M') 



(48) 



BARE VERTEX APPROXIMATION 



with EbvAo!,/3 given by Eq. (|47j). The inversion of 
Eq. jiq ) is given by Eq. (p9|), which is a set of four cou- 
pled integral equations for the four BVA Green functions, 
which must be solved simultaneously. 

From Eqs. (P2|) and (|23), the equations of motion for 
Xt(t) and the gap equation for x(t) is then given by: 



The bare vertex approximation (BVA) is obtained by 
setting the vertex function equal to its bare value in the 
exact equations for the one and two point functions. This 
is an energy conserving approximation which leads to in- 
tegral equations for the three- x vertex function as well as 
for the x-x-\ vertex function. The bare vertex approxi- 
mation consists of making the replacement 



r a ,i3, 1 [x](h,t2,t 3 ) = /q.,/3,7 Sc(ti,t 2 ) S c (ti,t 3 ) . (46) 



^+x(t)\x i (t) + ^K i (t,t)/i = 0, 



x(t) 



E 



1 



(49) 
(50) 



For the quantum roll, we further set Xi(t) = 0. This 
means that )Ci(t,t) = Ki{t,t) — 0, so that G a p(t,t') is 
diagonal, and results in the following set of equations for 
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FIG. 3: Vacuum graphs contributing to the 2PI part of the effective action ^[C?]. Solid lines represent the Gij(t,t') propagator, 
the wiggly to solid lines represent the ICi(t,t') and Ki(t,t') propagator, and wiggly lines are the T>(t,t') propagator. 



the Green functions: 
V(t, t') = D(t,t') 

- I dtif dt 2 £)(t,ti)n(ti,t 2 )©(t2 ! t / ) ! (5 1 ) 
Jc Jc 

g lJ (t,t') = G tJ (t,t') 

Yl I dt 4 d ^ G ik(t,ti)^ki(ti,t 2 )Qij(t 2 ,t') , 
,., Jc Jc 



(52) 



where 



x ij (t,t') = -g ij (t ) t')v(t',t). 



The gap equation for x{t) becomes: 

x(*) = §{^E^(m)/* 



(53) 



(54) 



In addition, for this case, the initial conditions imply that 
we can take Gij(t,t') and Gij(t,t') to be diagonal, which 



greatly simplify the integral equations. The BVA for the 
quantum roll requires that we solve equations (|5lj), (|52]), 
(p3|), and (^4j) simultaneously using the numerical meth- 
ods described in refs. |24| and p5fl . 

Because of the interest in using the BVA approxima- 
tion in QED (and QCD) plasma conductivity problems, 
we will discuss the integral equation one obtains for the 
vertex function in what follows. It was precisely because 
this approximation gives the sum of the graphs used in 
non-relativistic plasmas (see Fig. ^|) in conductivity cal- 
culations which gave both accurate results as well as giv- 
ing physical answers that initially interested us in this 
approximation. 

The three-point vertex functions for the BVA are given 
by functional differentiation of the inverse two point func- 
tions: 



A{h,t 2 , t 3 ) = r 000 (t u t 2 , t 3 ) = —„ ff 1 '* 2 ) (55) 



^ij(ti,t2,ts) = Fyo(il) t 2 , ts) 



5x(t 3 ) 



(56) 



and obtain the coupled integral equations: 



Hh,t 2 ,h) = -— / dt 4 / dt 5 G ik (ti , tj ) Yki {tj , h , h ) Qi j ft 5 , h ) Q ji (t 2 ,h) 



and 



(57) 



Yij{ti, t2, h) = 5ij S c (h,t 2 ) 6c(ti,ta) 

- J dt 4 J &h\^g ik {t l M)?ki{U,hM)Gi ] {h,t2)V{t 2 ,t l ) + g i] {t 1 M)^^^ (58) 

This is shown diagrammatically in Fig. ||. Looking at the diagrams, if we iterate these equations, we sum all the 
"rainbow" diagrams. As advertised, comparing these graphs with those shown in Fig. [IJ /Ci is approximated in the 
BVA by x exchange and K, 2 by x exchange in the i-channel. 

Let us show that this approximation is easy to obtain from the CJT formalism once we treat Q and T> and K, on 
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exactly the same footing. We choose for our approximation to T 2 [G] the 2-PI graphs shown in Fig. ||. This gives: 

A A ** C w C A A C •* C 



Since the T> propagator sums the contact term plus all the 
polarization bubbles II of the original quartic interaction 
gx 4 , if we reexpand I? in a power series in II then the 
first two terms in the series give the graphs used in the 
approximation ofjlj] andjl^]. The CJT action is given by 
Eq. (Ea). The stationary condition for G a a(t,t') gives: 



ST[x, G] 



sr 2 [G] 

SGaP 



= 0. 



or 



where: 

§r 2 [G] 



Sbvao,^[0](* ) * / ) = -2i 



5G a p(t,t>) 



(60) 



(59) 



r 



Carrying out the derivatives of T 2 [G] given in Eq. (p3S|), 
we find that ^BVAa,d{t,t') is exactly the same as found 



in Eq. (47) using the Schwinger-Dyson equations in the 
BVA approximation. The stationary condition for x a 
also gives the same equations of motion for Xj(t) and 
gap equation for \(t) as found in Eqs. (|4^) and (M) us- 
ing the Schwinger-Dyson equations in the BVA. Thus we 
conclude that the CJT action, as given in Eqs. ( (45| ) and 
(|59|) , gives exactly the same set of equations as in the 
Schwinger-Dyson BVA truncation. 



The energy for the BVA is obtained from 
if!) for the vertex function. We find: 



9h by using 



d 2 G(t, t')/i 



dtdt' 



+ x(t)G(t,t)/i-r 2 oX (t)-- 
t'=t 9 



X 2 (t) + j^V(t,t)/i 



[ dt i'D(h,t)Gi j (h,t)G j i(t,t 1 ) 



(61) 



where, for our case, we have set xi(t) = ii(t) = 0. Since 
the BVA equations are derived from an effective action, 
energy is conserved. 

VI. DYNAMICAL DEBYE SCREENING 
APPROXIMATION 

In plasma studies of the electric conductivity of fully 
ionized plasmas |2(], ^7j, it was found that in order to 
correctly determine the conductivity it was necessary to 
have an approximation where the photon propagator in- 
cluded the effects of dynamical Debye screening in the 
random phase approximation. This improved propaga- 
tor was then used in a scattering kernel in the kinetic 
equations. In our model, the \ field plays the roll of the 
photon in the dynamics of the Xi oscillators. The low- 
est approximation that includes the polarization effects 
in T> is precisely the leading order in large- TV approxi- 
mation to T>, namely T>q (see Eq. |6t]) which is discussed 
below in our derivation of the NLOLN approximation . 



The leading order in large- N approximation is similar in 
spirit to the random phase approximation. The equation 
for 1 (i, t') in leading order in large- A is given by: 

2V 1 (t, t') = -U c (t, + n (t, t') , (62) 

where 

+ J2 Xi(t) G l3 {t 1 t')x ] {t') . 

i,3 

In the QED plasma problem, the \ propagator be- 
comes the photon propagator and the delta function in 
T>o is replaced by the bare photon propagator. It is the 
bubble in n that leads to the Debye screening of the 
photon. It is because of our interest in QED that we call 
this approximation the DDSA. 
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Let us now specialize to the case when Xi(t) = 0. The The gap equation is: 
equation for the full x propagator Q is: 



g ij (t,t') = G ij (t,t') 

-V /"d*i/d*2Gifc(t,ti)E w (ti,t a )ftj(<a,0» ( 63 ) 
k,i Jc Jc 

with the self energy depending on the full Q and the lead- 
ing order in 1/N approximation to T> given by Eq. (|62|): 



z kl (t,t') = -g kl (t,t')D(t,t') 



(64) 



(65) 



There is a nontrivial vertex function in this approxima- 
tion given by: 



/C JC 

dt 4 / dt 5 K{U,t 5 ,t z ) D^^g^h) D{t 2 ,t b ) . 
c Jc 



(66) 



This equation can be obtained from the exact integral 
equation for T shown pictorially in Fig. ^ by making two 
approximations. The first is to approximate the exact 
three- x vertex function by the triangle graph, which is 
the leading term in the 1/N expansion of this function. 
The second is to replace the scattering kernels, Ki and 
K2 by single particle exchange in the t-channel. The rea- 
son for our studying this approximation is that, the same 
approximation made in QED can be shown to be the 
lowest order resummation scheme that preserves Ward 
identities (|9|). 

The DDSA approximation can be derived from an 
effective action by modifying slightly the approach of 
Cornwall, Jackiw and Tomboulis (CJT) pQ] . The discus- 
sion that follows here is due to Emil Mottola and Luis 
Bettencourt jl9| . Thinking of the fields x and x as P ar t 
of an N + 1 component field, and considering the case 
that (x(t)) = where there is no mixed propagator, one 
can write a CJT like action for the generating functional 
of the twice Legendre transformed effective action as: 

r[ X ,g,V] = S c uss[x] + '-TrilniV- 1 }} 

+ ^Tr{ In [g- 1 ]} + ^T T {V^ V+G^x] G-1} + T 3 [Q\. 

(67) 



here G' 1 (t, t') is defined by @ and V (t, t') by Eq. @. 
T>o(t,t') is considered an external parameter, and is not 
varied to obtain the equations of motion. In the DDSA, 
the 2-PI contribution to the action, ^[£7], for the case 
when Xi(t) — 0, is given by Eq. (^9[) with V(t,t') set 



equal to its LOLN value T>o(t,t'): 



r ? \g] = 



~ Ia? J2 J dtl J dh v o(h,t2) Qa (ti, h) g Jt (t 2 , h) 



(68) 



By varying the action (|67J), we reproduce Eqs. ( |63[ ) and 
(|65|). Although there is an effective action for the DDSA 
approximation, since T> is treated as an external time- 
dependent propagator, energy conservation is violated at 
order 1/N. At modest N we will find that this causes 
this approximation to become inaccurate after several os- 
cillations. However, it is more accurate at these modest 
values of N than the LOLN approximation, as well avoid- 
ing the unboundedness of the NLOLN approximation we 
discuss next. 



VII. THE LARGE- iV APPROXIMATION 

The large- A*" expansion is obtained from Eq. (|l3|) by 
first integrating over all the X4 and then evaluating the 
remaining functional integral over \ by steepest descent. 
The effective action, as a power series in 1/N, is ob- 
tained from the first Legendre transform of the generat- 
ing functional. In a previous paper ||, we obtained equa- 
tions for the next to leading order large-A approxima- 
tion (NLOLN) to the action, and gave numerical results 
for the quantum roll. For completeness, we review those 
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equations here. To order l/N, we obtain: 

rLarge-N = ^classf^] 

dttt^\n[G«%t)] + ±\n[Vo%t)] 

C ^ i 

where SciassM is given by Eq. (^Tj), and 2?^" 1 (t,t') is the 
inverse propagators for \ in lowest order in the l/N ex- 
pansion, given by 



Z>o (*, t') = D- l {t, t') + n (i, t') , (69) 



( 2) 

and where the second order Xi propagator Q\- ' [t, t) and 
self energy T,ij(t,t') to order l/N is given by: 



V [ dh[ dt 2 G ik (t,t 1 )X kl (t 1 ,t 2 )G lj {t 2 ,t'), (73) 



fe.i 



where 



s y (*» = 77 G y (*> W, t') - Xi (t) D (t, »j (t') . 



with 



53 &<(*) Gij{t,t') Xj {t') 



(70) 



Here D~ 1 (t, t') and Gl^{t,t') are the same as Eqs. @ 
that we defined earlier. 

The equations of motion for the classical fields Xi(t), 
to next to leading order in l/N, are given by: 



di 2 



-H 



J dt'G ij (t,t')V (t,t')x j (t') = 0, (71) 



with the gap equation for \(t) given by 

^) = |{e(^(*)+^E^ ) (^)a)-^} 



(72) 



We see here that the equation for Q is the expansion of 
the BVA equation in a series of l/N, truncated at first 
order. 

Let us now specialize to the case of the quantum roll 
problem where Xi(t) = 0. In that case the two point 
inverse propagator for the x field is 

G « M(*i,* 2 )- 5xi{tl)5xj{t2) 



with 



G^[x]{t l ,t 2 ) + Y lij {x]{hM), 



S y [x](t,0 = -G l3 (t,t')V (t,t') 



However it is G^' (t, t') which enters into Eq. ( [72] ) and not 
Gij{t,t'). Thus the solution for Gij(t,t'), which we might 
interpret as (xi(t)xj(t')}, does not enter into the dynam- 
ics of the solution! This Gu(t,t) is positive definite, but 
quickly blows up. 

The vertex function Tij(t±, t 2 ,ta) is given by: 



= Sc(ti,t 2 )5c(t 2 ,t 3 )5ij - — Gij(fi,f 3 ) Gji(ts,t 2 )T> (t 2 ,ti) 

- jj J c dnJ^dt 5 G ij (t 1 ,t 2 )V (t 1 ,U)A {U,t 5 ,t3)Vo(t 5 ,t 2 ) , 



IC JC 

where the lowest order in l/N 3-x vertex is given by 

SV^(t 4 ,t 5 ) 



Ao(t4,*5,*3) 



= — Jr E ^J'C* 4 ' ^(*3) *5) Gli(t 5 , t 4 ) . 



We immediately see that this is not an integral equation but again, is the lowest order in l/N contribution to Eq. (|57|). 
The inverse x propagator gets 1/iV corrections which are of two types, one is a self energy correction to the x 
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propagator and the other is a new three loop graph containing two lowest order \ propagators. We find 

t^IYarge-N [x, x] 



V-\l,2) 



--5 c (l,2)-no(l,2)- V / dtJ dt i G ij {t l ,h)T, jk {t s ,U)Gu{t^h)G u {h,t l ) 
9 ijk i J c J c 

+ I At J dt 4 f dtJ dt 6 Ao(t 4 ,t 1 ,t 3 )Vo(t3,t 5 )A (t 5 ,t 2 ,t 6 )Vo(t 6 ,t 4 ) . 

J C J C J C J c 



r 



The last term in this equation is a 1 /N correction to the 
vertex function. However, it is T>q and not T> which en- 
ters Eq. (|73|), so that the BVA and the 1/N expansion 
will differ only by terms of order 1/N 2 . The BVA ap- 
proximation treats x and x on exactly the same footing, 
whereas the large- TV expansion treats x exactly, but then 
expands in loops of x- So at order 1/N 2 , the large- N ap- 
proximation will contain graphs omitted from the BVA 
approximation, and vice-versa. 



VIII. RESULTS AND CONCLUSIONS 

In this section we present the results of exact numer- 
ical simulations of the quantum roll, using initial con- 
ditions described in our previous paper on the large-V 
approximation We choose as our dimensional mass 
scale the second derivative of U(r) at the minimum of 
the effective one dimensional potential U(r). This mass 
scale was chosen to have value M 2 = 2. In terms of this 
mass scale, the coupling constant as well as the rescaled 
ro are of order one for all N. The exact manner in which 
g and ro runs with N is described in ref . . 

As N — > oo the Hartree and leading order large- TV ap- 
proximation become exact and an initially gaussian wave 
packet remains gaussian with width equal to (x 2 (t)} os- 
cillating in a known manner. At modest N, 10 < N < 20 
an initially Gaussian wave function develops a large num- 
ber of nodes and so the wave function even at modest 
times is of the form Gaussian time a high order polyno- 
mial. In spite of this, (x 2 (t)) shows rather simple be- 
havior. It oscillates with a constant amplitude for a rea- 
sonable period of time with an envelope that oscillates 
with a much longer time constant which increases with 
N. The Hartree and leading order large- N approxima- 
tions just oscillate with fixed amplitude. The NLOLN 
blows up in this regime. BVA attempts to track the con- 
traction of the envelope but then contracts to a fixed 
point. The DDSA violates energy conservation at order 
1 /N so it becomes numerically inaccurate when 1 /N ef- 
fects become important which is at a time t ex. N. Both 
BVA and DDSA do however stay bounded and positive 
definite during the time period of our numerical simu- 
lations. Higher order correlation functions show more 
complicated behavior and the approximations presented 
here are only accurate for a few oscillations in the regime 



3 < N < 20 consistent with the increasingly complicated 
evolving structure of the wave function. 

In Figs. H through^, we show the results for (x 2 (t)) as 
a function of t, comparing the bare vertex, the dynamic 
Debye screening, and the large- ./V approximations to the 
exact solution, for N = 3, 10, and 21. In Figs. ^ to ||, we 
show the same results for (x(t)) as a function of t, and 
in Figs. ^ through [Tl|, we give the results for (x 2 (t) [For 
detailed views of these figures in color, see our web site 



at : http : / / www. theory. unh . edu/resum | . 

In our previous studies || of the large- N approxima- 
tion, we found that the next to leading order large- N ap- 
proximation had the feature that the effective potential 
was not defined at small x for N < 20, for our param- 
eter set, and it was not until N was greater than about 
20 that the large-V expansion produced bounded values 
for (x 2 (t)). This result is reproduced here. For the limit 
N — > oo the quantity (x 2 (t)) corresponds to harmonic 
oscillations. At finite N, however, the exact solution for 
N > 21 has the property that the envelope of these oscil- 
lations contracts. As noted in the figures, only the bare 
vertex approximation attempts to follow this contraction. 
At N = 21, the BVA is accurate up to a t « 130 before 
overshooting and then oscillating about a fixed point. 
This fixed point behavior shows that this approximation 
still neglects some important quantum phase information 
present in the exact solution. 

In contrast to the NLOLN approximation, which 
breaks down for V < 21, both the BVA and the DDSA 
have the feature that (x 2 (t)) remains positive definite, as 
well as being bounded at all TV. This is true for all the 
expectation values that contribute to the energy. This 
conclusion is purely based on numerical evidence. We do 
not have a proof that this approximation corresponds to 
a positive definite probability distribution. However, all 
the moments we have studied (a total of five, as shown 
in Fig. |l2|), are all bounded. 

The DDSA is more accurate than the second order 
large-V approximation for N less that 20, but for N 
greater than 20, the reverse becomes true. However, nei- 
ther approximation captures the true nonlinear shrinking 
of the envelope of the oscillations, even for N greater than 
20. 

Energy is conserved for the bare vertex and the second 
order large-V approximations, but not for the dynamic 
Debye screening approximations, as pointed out in sec- 
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N = 3 



exact 
Bare Vertex 
Dynamic Debye Screening 
2nd order large-N 
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t 

FIG. 4: Plot of (x 2 (t)) as a function of t, comparing the bare vertex, the dynamic Debye screening, and the large- N approxi- 
mations to the exact solution, for N = 3. 




FIG. 5: Plot of (x 2 (t)) as a function of t, comparing the bare vertex, the dynamic Debye screening, and the large- N approxi- 
mations to the exact solution, for N = 10. 
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FIG. 6: Plot of (x 2 (t)) as a function of t, comparing the bare vertex, the dynamic Debye screening, and the large- N approxi- 
mations to the exact solution, for N = 21. 
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FIG. 7: Plot of (x{t)) as a function of t, comparing the bare vertex, the dynamic Debye screening, and the large-iV approxi- 
mations to the exact solution, for N = 3. 
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FIG. 8: Plot of (x{t)) as a function of t, comparing the bare vertex, the dynamic Debye screening, and the \arge-N approxi- 
mations to the exact solution, for N = 21. 
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FIG. 9: Plot of {x 2 (t)} as a function of t, comparing the bare vertex, the dynamic Debye screening, and the large- ./V approxi- 
mations to the exact solution for N = 3. 
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N = 10 
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FIG. 10: Plot of {x 2 {t)) as a function of t, comparing the bare vertex, the dynamic Debye screening, and the large- iV approxi- 
mations to the exact solution for N = 10. 
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FIG. 11: Plot of { X 2 (t)) as a function of t, comparing the bare vertex approximation to the exact solution for N = 21. 
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FIG. 12: Plot ol various contributions to the energy for the bare vertex approximation as a lunction of t for TV = 10. 



tion VI . This is a serious drawback to the dynamic Debye 
screening approximation. 

In all these figures, one can see that the bare vertex 
approximation tries to follow the envelope of the exact 
curve, whereas the dynamic Debye screening approxima- 
tion does not do so. This is particularly striking for the 
cases when TV is less than 21, where the dynamic Debye 
screening approximation yield unphysically large values 
for the expectation values. 

In the BVA approximation we observe that (x 2 {t)) at 
late times has an envelope of decreasing oscillations about 
a fixed point. In fact as seen in Fig. [lj all the contribu- 
tions to the energy in the BVA have the same feature 
that they asymptote to a fixed point. In Fig. [ij we dis- 
play all five contributions to the energy at N = 10 to 
demonstrate this fact. In contrast, as seen in the very 
long time run shown in Fig. [ll], the exact solutions ex- 
hibit "recurrence" patterns of motion which are not cap- 
tured in the BVA. In the 1 + 1 dimensional field the- 
ory simulations of ref . jl5| , all the Fourier components of 
the two particle correlation function showed this behav- 
ior which was given as evidence for thermalization. So 
one hopes that this "defect" of the BVA approximation 
in a quantum mechanics setting, will instead have the 
correct physics of thermalization in a field theory appli- 
cation where Poincare recurrence times are expected to 
become very large. To see if this is true, we intend to 
study the BVA in classical 1+1 dimensional field theory 
where again exact simulations can be performed ||. 

In summary we have found that both resummation 



methods described here, the BVA and the DDSA, pro- 
duce positive definite and apparently bounded results for 
expectation values at all values of N. The bare vertex 
approximation appears to provide the best description of 
the motion, but cannot describe recurrences of the mo- 
tion. Still, it provides an energy conserving and reason- 
ably accurate description, and is a dramatic improvement 
over the next to leading order large-iV approximation 
when TV < iV cr it = 21. As mentioned earlier, in the single 
particle quantum mechanics problem we studied here, the 
graphs do not correspond to particle collisions, so there 
is no possibility of studying thermalization. Thermaliza- 
tion questions need to be addressed in field theory appli- 
cations. It will be important to show that the BVA ap- 
proximation will lead to thermalization of arbitrary ini- 
tial data as found in the 3-loop approximation of ref. |Ifi[] 
when applied to 1+1 dimensional quantum field theory. 
We would also like to study the analogue of the BVA 
approximation for a gaussian ensemble of initial condi- 
tions for a 1+1 dimensional classical field theory since 
that can also be studied exactly numerically ||. These 
authors have shown that the classical field theory indeed 
thcrmalizes and we would like to know how accurately 
the classical version of our approximation captures this 
physics. This will be the subject of a future publication. 
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